ms  Drru' 
b*  copy  ?u  ,.i  - 

IGMIFICANI  NEVCKR  OF  FAG 
RODUCE LEGIBLE  „> 


The  Florida  State  University 

Department 
of 

Statistics 

Tallahassee.  Florida 


X I  CARLE. . 

AiTAINSD  A  / 

LiYHJWHCOgOI  ' 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


D  D  G 

iElSEIlflilE 


NOV  7  W® 

tei^binrEi^ 

E 


A  MULTIVARIATE  METHODOLOGY  FOR  THE  ANALYSIS  OF 
WEATHER  MODIFICATION  EXPERIMEI; 


/O 


by 

Elton  Scott 


FSU-Statistics  Repert-No.  M514 
~ONR  Technical  Report  No.  147' 


/ 


>{ 

>/ 

AugO— ^u'lD79 

The  Florida  State  University 
Department  of  Statistics  / 
Tallahassee,  Florida  32306 


Thia_woxk_was  supported  by  the  Office  of  Naval  Research  Contract  No. 
N00014-76-C-0394k/  Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government. 


PREFACE 


The  main  body  of  this  technical  report  is  the  paper  "A 
Multivariate  Methodology  for  the  Analysis  of  Weather  Modification 
Experiments".  This  paper  will  be  printed  in  the  Proceedings  of 
the  Workshop  on  Statistical  Design  and  Analysis  of  Weather  Modifi¬ 
cation  Experiments,  Tallahassee,  October  1978.  In  addition  to 
that  paper  this  report  includes  an  Appendix  that  provides  supple¬ 
mental  information  on  the  principal  components  results  summar¬ 
ized  in  the  main  body.  The  Appendix  gives  details  on  station  lo¬ 
cations,  estimated  station  precipitation  means  and  standard  devi¬ 
ations,  eigenstructure  estimates,  as  well  as  station-to-station 
and  station-to-component  correlations. 

A  complete  list  of  technical  reports  on  this  contract  is 
appended . 
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ABSTRACT 

This  paper  develops  applications  of  multivariate  statistical 
models,  particularly  principal  component  analysis,  to  the  analysis 
of  data  from  weather  modification  experiments.  The  efficacy  of 
these  multivariate  applications  is  examined  by  applying  the  proposed 
models  to  data  from  Phase  I  (1967-71)  of  the  Santa  Barbara  Convec¬ 
tive  Band  Seeding  Program  conducted  for  the  Navy  by  North  American 
Weather  Consultants.  Multivariate  summary  measures  of  precipitation 
are  developed  and  multivariate  methods  are  given  to  analyze  the  effects 
of  cloud-seeding  on  precipitation.  Results  from  these  models,  based 
on  the  above-mentioned  data  set,  are  reported  along  with  conclusions  and 
suggestions  for  further  work.  An  Appendix  provides  detailed  summary 
statistics  for  the  analyses. 


*This  research  was  supported  by  the  Office  of  Naval  Research  under 
Contract  No.  N00014-76-C-0394.  Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of  the  United  States  Government. 


1.  INTRODUCTION 


Weather  modification  experiments  generally  produce  multi¬ 
variate  data.  Both  the  precipitation  measurements  and  concomi¬ 
tant  variables  for  a  given  experimental  unit  are  usually  repre¬ 
sented  as  a  vector  of  measurements,  and  as  such,  multivariate 
methods  are  appropriate  for  the  analysis  of  the  data. 

Logically,  the  data  analysis  can  proceed  in  two  steps. 

First,  the  precipitation  measurements  must  be  represented  by  an 
appropriate  summary  measure  and,  second,  these  summary  measures 
for  the  experimental  units  can  be  used  to  examine  the  effects 
of  seeding. 

Below,  in  Sections  (2)  and  (3),  models  for  these  two  stages 
are  reviewed  and,  at  each  stage,  the  models  are  applied  to  the 
Phase  I  data  from  the  Santa  Barbara  Convective  Seeding  Program 
(SBA-I  data  hereafter)  (Thompson,  Brown,  and  Elliott,  1975). 

The  final  section  of  the  paper  gives  conclusions  on  the  analyses, 
limitations  of  this  application,  and  suggestions  for  further 
work. 

2.  SUMMARY  MEASURES  OF  PRECIPITATION 

With  most  weather  modification  experiments,  a  large  net¬ 
work  of  raingages  is  set  up  over  designated  areas  to  collect 
precipitation  measurements.  Often  measurements  from  each  ex¬ 
perimental  unit  are  simply  averaged  and  the  simple  average  is 
taken  to  represent  the  volume  of  precipitation  on  the  area  for 
that  experimental  unit.  Inasmuch  as  the  raingages  are  irreg¬ 
ularly  spaced,  the  measurements  include  some  noise,  and  the 
measurements  tend  to  be  positively  (but  imperfectly)  correlated, 
the  adequacy  of  such  averages  to  represent  the  full  data  set 
can  be  questioned. 

This  section  provides  a  brief  intuitive  description  of  a 
multivariate  statistical  model  for  dealing  with  correlated  data 
sets.  Principal  Components  Analysis  (PCA  hereafter).  A  method 
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that  uses  PCA  to  produce  multivariate  summary  measures  of  pre¬ 
cipitation  data  from  the  experimental  units  is  presented  and 
other  uses  of  PCA  are  considered.  This  section  close  with  an 
evaluation  of  results  obtained  when  these  methods  were  applied 
to  the  SBA-I  data.  The  measures  produced  at  the  first  stage  can 
be  used  in  the  second  stage  to  examine  the  effects  of  seeding 
and  the  second  stage  is  considered  in  Section  (3). 

2.1  The  Concept  of  PCA 

The  data  array  is  arranged  so  that  the  vectors  of  precipita¬ 
tion  measurements  for  the  experimental  units  (convective  bands) 
are  regarded  as  the  observation  vectors  and  the  measurements  at 
raingage  stations  as  (possibly)  correlated  variables.  Let  y. 
be  the  measurement  at  Station  i  in  the  ot  experimental  unit. 

We  use 


Zia=  (yia  "  yi)/si  1  =  P5  3  =  1,2,. ..,N, 

N 

to  calculate  the  p  x  p  correlation  matrix,  R  =  [  7  z.  z.  /N] . 

r  r  ~  L  **,  ia  ia J ' 

o=l  J 

where  y^  is  the  mean  over  the  N  observations  at  Station  i,  and 
s^  is  the  corresponding  estimated  standard  deviation  calculated 
from  the  N  values  (observed  at  Station  i  for  the  N  convective 
bands) . 

The  principal  components  to  be  determined  depend  on  the 
characteristic  roots  and  vectors  of  R.  Let  the  ktl1  character- 
istic  root  and  corresponding  characteristic  vector  be  given  by 
X^.  and  where  elements  of  ^  are  a^,  i  =  l,2,...,p,  k  =  1, 

2,...,p.  The  k**1  principal  component  for  the  <xt^  experimental 
unit  is: 


Cka  ~.^.aki  Zict‘ 
i=l 


Of  course  the  ^  and  X^  must  be  estimated  and  these  estimates 
(designated  as  ^  and  X^  below)  are  used  to  estimate  the  principal 


component  values  for  each  observation  vector. 

The  estimated  product -moment  correlation  between  the  i^ 
variable  and  the  kth  component  is  given  by  ^i1^^  (Morrison, 

1976  (p.  271)).  Thus,  coefficients  estimated  from  standardized 
data  (as  above)  indicate  the  extent  and  sign  of  the  association 
between  variables  (stations)  and  principal  components. 

The  principal  components  are  ordered  by  the  proportions  of 
total  variance  represented  by  them.  If  the  characteristic  roots 
are  extracted  from  a  correlation  matrix,  as  in  the  present  appli¬ 
cation,  the  sum  of  the  characteristic  roots  will  be  p  and  the 
proportion  of  the  total  variance  attributable  to  the  kl  component 
will  be  A^/P- 

Since  PCA  is  applied  so  that  the  dimensionality  of  response 
vectors  can  be  reduced,  one  would  like  the  first  few  principal 
components  to  account  for  most  of  the  variance  of  the  data.  Thus, 
if  J  of  the  components  "adequately"  represent  the  data,  the  p- 
dimensional  variate  vector  could  be  summarized  by  the  corres¬ 
ponding  J-dimensional  principal  component  vector. 

All  variables  usually  have  high  pairwise  correlations  with 
the  first  principal  component  and  further,  the  second  and  subse¬ 
quent  components  may  exhibit  strong  associations  with  subsets 
of  variables.  If  the  first  few  principal  components  account 
for  most  of  the  variance,  and  if  each  principal  component  is 
highly  correlated  with  one  or  more  representative  stations,  future 
data  collection  might  be  simplified.  Data  from  these  representa¬ 
tive  stations  should  essentially  contain  the  information  of  the 
principal  components  with  which  they  are  associated,  and,  since 
the  components  represent  the  essential  information  from  the 
full  network  of  stations,  most  of  the  information  from  the  full 
network  should  be  reflected  in  the  information  from  the  represent¬ 
ative  subset.  If  these  assumptions  hold,  the  efficiency  of  future 
data  collection  could  be  improved  since  costs  should  be  reduced 
by  improvement  in  the  reliability  of  data  collection  at  the 
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subset  of  stations  and  collection  of  data  from  the  representative 
subset  of  stations  only,  rather  than  from  the  full  network. 

2.2  Application  of  PCA  to  Weather  Modification  Data 

The  application  of  PCA  to  weather  modification  experiments 
would  proceed  with  initial  estimation  of  principal  components 
from  available  data.  Given  that  the  estimates  indicate  that  the 
first  few  principal  components  do  account  for  most  of  the  vari¬ 
ance  in  the  full  data  set,  a  subset  of  key  stations  could  be 
selected  to  reduce  the  costs  of  future  data  collection  in  weather 
modification  experiments.  If  the  data  are  pooled  in  this  manner, 
we  must  assume  that  the  principal  axes  of  data  are  not  affected 
by  seeding.  That  is,  the  approach  assumes  that  the  principal 
components  are  the  same  for  both  seeded  and  unseeded  experimental 
units. 

PCA  can  also  be  applied  to  evaluate  weather  modification  ex¬ 
periments.  The  estimates  of  principal  component  values  for  exper¬ 
imental  units  are  vector  summary  response  measures.  These  vectors 
can  be  compared  to  other  summary  measures,  such  as  the  mean  pre¬ 
cipitation,  to  determine  the  effectiveness  of  the  PCA  approach 
to  summarization  of  the  data.  If  this  approach  is  useful,  the 
vectors  for  the  seeded  data  could  be  compared  to  the  vectors 
from  the  unseeded  data  to  determine  the  effects  of  seeding.  This 
analysis  is  pursued  in  Section  (3)  of  this  paper.  The  patterns 
of  the  station-component  correlations  may  also  prove  useful  in 
evaluation  of  the  nature  of  the  precipitation  response  over  an 
area. 

In  most  applications  of  PCA,  the  simple  mean  of  the  variables 
is  associated  with  the  first  principal  component.  If  the  band 
means  and  the  first  principal  component  scores  are  highly  corre¬ 
lated,  then  the  second  and  subsequent  components  must  repre¬ 
sent  essential  features  of  the  variates  that  are  not  reflected 
in  their  mean.  For  precipitation  data,  the  second  and  subsequent 
components  could  be  associated  with  subarea  concentrations  of 
precipitation. 
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2.3  Evaluation  of  PCA  Results 


An  obvious  question  on  the  above  procedure  is  "How  well  do 
the  first  few  components  represent  the  data  from  the  full  network 
of  stations?"  This  question  can  be  answered  through  use  of  re¬ 
gression  analysis  to  "predict"  summary  response  measures  of  inter¬ 
est  with  the  principal  component  values  as  the  independent  vari¬ 
ables.  Thus,  if  the  mean  is  the  summary  measure  of  interest,  it 
should  be  reasonably  represented  as  a  linear  function  of  the  first 
few  principal  components. 

An  alternative  evaluation  of  the  usefulness  of  the  PCA  ap¬ 
proach  depends  on  the  ease  of  interpretation  of  the  components. 

The  station-component  correlations  may  be  interpretable  through 
examination  of  a  graphical  presentation  of  the  station-component 
correlations.  Interpretation  procedures  often  are  tenuous  how¬ 
ever,  because  the  principal  components  sometimes  do  not  present 
interpretable  patterns  for  the  variables. 

If  the  PCA  results  compare  favorably  to  other  summary  meas¬ 
ures,  the  usefulness  of  a  key  subset  of  stations  in  representing 
the  summary  measure  can  also  be  evaluated  from  the  regression 
of  the  summary  measure  on  the  key  station  measurements. 

Thus,  the  steps  suggested  in  application  of  PCA  to  precipi¬ 
tation  data  are: 

1.  Estimate  principal  components  over  the  full  network  of 
stations. 

2.  Examine  the  efficacy  of  the  principal  components  as  a 
summary  response  measure. 

3.  Interpret  graphical  presentations  of  the  station- compon¬ 
ent  correlations  as  appropriate. 

Further  steps  to  reduce  future  data  collection  costs  would  be: 

4.  Select  a  key  subset  of  stations  that  are  highly  corre¬ 
lated  with  the  first  few  principal  components. 

5.  Test  the  effectiveness  of  these  key  stations  in  repre¬ 
sentation  of  the  information  in  the  summary  measure  of  interest. 
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Final  steps  would  then  use  the  above  results  to  test  for  the 
effects  of  seeding.  These  testing  steps  are  detailed  in  Section 
(3)  of  this  paper.  The  remaining  subsections  of  this  section 
report  the  results  of  the  above  steps  for  data  from  Phase  I  of 
the  Santa  Barbara  experiments  (Thompson,  Brown,  and  Elliott; 

1975). 

2.4  The  Santa  Barbara  Data 

The  Santa  Barbara  experiments  and  data  are  briefly  reviewed 
by  Bradley,  Srivastava,  and  Lanzdorf  (1977)  and  detailed  in 
Thompson,  Brown,  and  Elliott  (1975).  The  frequency  of  missing 
data  at  many  of  the  stations  and  the  nature  of  PCA  required  an 
adjustment  in  the  set  of  stations  used  in  this  analysis.  The 
basic  areas  of  interest  are  the  "Control"  and  "Target"  areas 
(see  Bradley,  Srivastava,  and  Lanzdorf,  1977),  but  the  number  of 
stations  in  each  area  was  reduced  considerably  because  of  missing 
data. 

The  control  and  target  areas  were  analyzed  separately  to  be 
consistent  with  the  response  surface  approach  of  Bradley,  Srivas¬ 
tava,  and  Lanzdorf  (1977).  The  resulting  principal  component 
values  could  be  used  as  covariates  or  summary  measures  in  subse¬ 
quent  analyses. 

As  mentioned  above,  many  stations  had  large  numbers  of  mis¬ 
sing  measurements  (a  complete  set  of  measurements  for  a  station 
would  include  precipitation  measurements  for  all  107  Phase  I  bands). 
Stations  with  large  numbers  of  missing  values  present  the  possi¬ 
bility  of  a  sample  covariance  matrix  that  is  not  positive  defi¬ 
nite.  When  PCA  was  applied  to  all  stations,  irrespective  of  mis¬ 
sing  data,  sizeable  negative  characteristic  roots  occurred,  es¬ 
pecially  for  the  target  area.  To  avoid  this  problem,  all  stations 
with  more  than  22  missing  band  values  were  exluded  from  the  analy¬ 
sis  (this  cutoff  point  is  arbitrary,  but  it  did  improve  the  be¬ 
havior  of  the  latent  roots).  This  deletion  procedure  yielded 
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nineteen  control  stations  (fifteen  were  deleted)  and  seventy-two 
target  area  stations  (thirty-eight  were  deleted).  A  complete  list 
of  included  stations  is  given  in  the  Appendix. 

The  elimination  of  stations  could  have  systematically  ex¬ 
cluded  some  geographic  subareas  in  either  Target  or  Control  Areas. 
To  test  for  this  possibility,  multivariate  analysis  of  variance 
(MANOVA)  was  applied  to  the  vector  of  location  coordinates  (lat¬ 
itude,  longitude,  and  altitude).  MANOVA  is  the  multivariate  ana¬ 
log  of  univariate  ANOVA  (Morrison,  1976).  The  "treatment"  groups 
were  the  included  and  excluded  sets  of  stations,  and  the  null  hy¬ 
pothesis  was  that  there  were  no  differences  in  the  mean  location 
vectors  foT  the  two  groups. 

The  null  hypothesis  could  not  be  rejected  for  the  Control - 
Area  groups  but  for  the  Target- Area  groups,  the  null  was  rejected 
(a  =  .05).  This  implies  that  there  are  significant  differences 
in  the  coordinates  of  the  excluded  and  included  stations.  This 
could  produce  systematic  variations  for  results  on  the  included 
subset  relative  to  results  for  the  full  network  of  stations. 

The  excluded  stations  for  the  Target  Area  were  primarily  in 
the  northeast  quadrant  of  the  Target  Area,  so  they  are  remote 
relative  to  the  seeding  site.  Because  these  stations  are  gener¬ 
ally  remote,  the  systematic  exclusion  of  these  stations  should 
not  limit  the  extent  to  which  results  below  can  be  generalized. 

In  any  event,  the  advantages  of  including  stations  with  large 
amounts  of  missing  data  were  not  considered  to  be  important 
enough  to  overcome  the  problems  produced  by  including  these  sta¬ 
tions. 

2.5  Principal  Component  Estimates 

As  specified  above,  the  data  were  standardized  and  correla¬ 
tion  matrices  were  estimated  for  the  stations  over  the  107  bands. 
The  pairwise  simple  correlations  estimates  on  R  were  based  on  all 
data  available  in  common  for  each  pair  of  stations.  Character¬ 
istic  roots  and  associated  characteristic  vectors  were  obtained 
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separately  for  the  72-square  correlation  matrix  for  the  Target 
Area  and  for  the  19-square  correlation  matrix  for  the  Control 
Area  and  corresponding  values  for  the  principal  components  were 
obtained  for  each  convective  band.  Table  1  gives  the  character¬ 
istic  roots,  proportions  of  variance,  and  cumulative  proportions 
of  variance  that  are  accounted  for  by  the  first  few  principal  com¬ 
ponents  of  both  Target  and  Control  Areas.  These  results  suggest 
that  the  first  few  principal  components  do  indeed  account  for 
much  of  the  variance.  Thus,  a  few  dimensions  of  the  principal 
component  vector  accounts  for  most  of  the  variance  in  the  precip¬ 
itation  measurement  vectors. 

TABLE  1:  EIGENVALUES  AND  PERCENT  VARIANCE  ASSOCIATED 
WITH  PRIMARY  PRINCIPAL  COMPONENTS 


Princ.  Component 
Number 

Eigenvalue 

%  of  Var. 

Cum.  %  Var. 

Target  Area 
(72  Stations): 

1 

51.33 

71.3 

71.3 

2 

4.80 

6.7 

78.0 

3 

4.23 

5.9 

83.8 

4 

2.07 

2.9 

86.7 

5 

1.39 

1.9 

88.7 

Control  Area 
(19  Stations) : 

1 

14.45 

76.1 

76.1 

2 

1.27 

6.7 

82.8 

3 

.89 

4.7 

87.4 

4 

.59 

3.1 

90.5 

5 

.37 

1.9 

92.4 

2.6  Response  Information  from  Principal  Components 

The  information  content  of  the  first  few  principal  components 
of  the  data  is  analyzed  below.  Separate  band  precipitation  means 
(and  band  precipitation  standard  deviations)  were  estimated  from 
the  stations  included  in  the  PCA  data  over  Target  and  Control 
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Areas.  These  values  were  used  as  dependent  variables  and  were  re¬ 
gressed  on  the  first  few  principal  component  values  for  the  Tar¬ 
get  and  Control  Areas  (The  first  six  principal  components  for  the 
Target  Area  and  the  first  five  principal  components  for  the  Con¬ 
trol  Area).  Generally  means  and  standard  deviations  of  precipi¬ 
tation  measurements  for  experimental  units  are  positively  (but 
not  perfectly)  correlated.  It  is  possible  though,  for  cloud  seed¬ 
ing  to  affect  the  scale  parameter  of  the  distribution  without 
affecting  the  location  parameter,  so  standard  deviations  were  al¬ 
so  included  as  dependent  variables  in  similar  regressions.  The 
linear  regression  results  are  reported  in  Table  2. 

The  regression  results  confirm  that  the  first  few  components 
yield  near-perfect  predictions  of  the  means.  Although  it  is  not 
reported  in  Table  2,  the  first  principal  component  alone  accounted 
for  most  of  the  variance  in  band  means.  The  R  for  simple  regres¬ 
sions  of  the  means  on  the  first  principal  component  were  .997  for 
the  Target  Area  and  .985  for  the  Control  Area.  Thus  the  first 
principal  component  generally  reflects  the  summary  information 
in  the  band  means  and  the  second  and  other  principal  components 
reflect  further  summary  response  information. 

Values  of  R2  for  the  standard  deviations  are  lower  but  the  asso¬ 
ciation  with  the  first  few  principal  components  is  still  quite 
strong.  Here,  the  first  principal  component  is  again  strongly 
associated  with  the  standard  deviations  but  other  components  are 
also  important. 

The  statistics  given  in  Table  2  demonstrate  the  efficacy  of 
the  principal  components  as  summary  response  measures  for  precip¬ 
itation.  Response-surface  integrated  volumes  were  not  directly 
considered  because  these  values  are  highly  correlated  with  the 
simple  band  means  (Bradley,  Srivastava,  and  Lanzdorf;  1977). 

2,7  Interpreting  of  the  Principal  Components 

As  mentioned  above,  principal  components  are  not  always  in¬ 
terpretable  as  recognizable  constructs.  In  the  present  case 
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TABLE  2:  REGRESSION  RESULTS  PREDICTING  MEANS  AND 


T 


however,  the  results  do  yield  readily  interpretable  results. 

The  station-component  correlation  estimates  (r^  =  )  for 

the  first  three  components  of  both  Target  and  Control  Areas  are 
given  in  Figures  1  through  6.  The  left  side  of  the  horizontal 
bars  represent  station  locations  in  degrees  and  hundredths  and 
the  vertical  spikes  represent  the  r,  .  or  the  simple  correlation 
between  the  k  component  and  the  itfl  station  (the  scale  is  2 
inches  equals  a  correlation  of  1.0).  Correlations  between  sta¬ 
tion  precipitations  and  the  first  component  (Figures  1  and  4)  are 
large  for  all  stations  in  both  areas.  For  the  second  component 
(Figures  2  and  5),  the  correlations  are  associated  with  the  loca¬ 
tion  of  the  stations  relative  to  the  coast.  The  signs  are  oppo¬ 
site  for  the  inland  stations  relative  to  the  coastal  stations. 
Given  the  topography  of  the  Santa  Barbara  area  (i.e.,  mountain 
ridges  along  the  coast) ,  this  component  is  interpreted  as  an 
orographic  component  of  precipitation.  This  pattern  holds  for 
both  Target  and  Control  Areas,  as  does  the  third  component  pat¬ 
tern  (Figures  3  and  6) . 

The  magnitudes  and  signs  of  the  correlations  for  the  third 
component  vary  directly  with  the  East -West  location  of  the  sta¬ 
tions.  Again,  the  pattern  is  clearly  present  for  both  Target 
and  Control  areas.  The  patterns  for  the  other  principal  compo¬ 
nents  are  not  clearly  interpretable  and  are  not  displayed,  but 
these  components  are  relatively  unimportant  since  most  of  the 
variance  is  associated  with  the  first  three  components. 

2.8  Selecting  Key  Stations 

The  final  part  of  this  section  reports  on  the  use  of  PCA 
to  select  a  small  subset  of  key  stations  that  can  be  used  repro¬ 
duce  summary  response  measures. 

A  subset  of  stations  that  were  highly  correlated  with  the 
second  through  the  sixth  components  were  selected  (The  first 
component  was  not  considered  because  it  is  highly  correlated  with 
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FIGURE  2:  CORRELATIONS  BETWEEN  PRECIPITATION  AND  PRINCIPAL 

COMPONENT  2  BY  STATIONS,  TARGET  AREA. 
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FIGURE  3:  CORRELATIONS  BETWEEN  PRECIPITATION  AND  PRINCIPAL 

COMPONENT  3  BY  STATIONS,  TARGET  AREA. 
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most  stations).  For  the  Target  Area,  two  stations  were  selected 
for  each  of  the  five  components  and  for  the  Control  Area  only- 
one  station  was  selected  for  each  component.  Two  other  objectives 
were  considered  in  addition  to  maximization  of  the  correlation 
for  the  Target  Area  stations:  (i)  Stations  with  fewer  missing 
values  were  chosen  where  several  stations  had  nearly-equal  correl¬ 
ations;  (ii)  For  the  second  station  in  the  Target  Area,  the  sta¬ 
tion  was  also  chosen  to  be  geographically  separated  from  the 
first  station.  The  selected  stations  are  given  in  Table  3  along 
with  station-component  correlations.  In  both  the  Target  and  Con¬ 
trol  stations,  some  stations  were  highly  correlated  with  two  of 
the  components.  This  resulted  in  only  eight  key  Target  Area  sta¬ 
tions  and  four  key  Control  Area  stations. 

The  station-component  correlations  are  not  the  primary  con¬ 
cern  here.  These  key  station  measurements  should  reproduce  the 
summary  response  measures  if  this  approach  is  to  be  useful.  The 
results  for  a  series  of  regressions  on  these  key  stations  are 
reported  in  Table  3. 

For  the  Target  regressions,  the  values  of  the  simple  correla¬ 
tion  coefficients  are  good  (but  not  outstanding)  for  the  means, 
standard  deviations,  and  the  first  two  principal  components. 

This  suggests  that,  if  means  are  chosen  as  the  summary  measure, 
this  approach  may  be  useful  in  reducing  the  number  of  data  collec¬ 
tion  points  for  large  areas.  The  results  on  the  standard  devia¬ 
tion  and  the  second  through  the  sixth  component  are  not  adequate 
for  reliable  summarization. 

It  is  almost  certain  that  better  results  could  be  obtained 
here  by  the  inclusion  of  stations  that  are  highly  correlated  with 
the  first  principal  component  (several  station/first-component 
correlations  exceeded  .95).  All  stations  in  the  regressions  neces¬ 
sarily  had  lower-than-average  correlations  with  the  first  princi¬ 
pal  component  because  of  their  higher-than-average  correlations 
with  components  orthogonal  to  the  first  component.  Again,  this 
conjecture  was  not  tested  here  but  should  be  considered  in  any 


19 


TABLE  3:  STATIONS  SELECTED  TO 

REPRESENT  PRINCIPAL  COMPONENTS 


TABLE  4:  REGRESSIONS  ON  KEY  STATIONS 


Sum. 

Stat.  and 

Ind. 

Var.  Mean 

Std.  Dev. 

PCS1* 

PCS2* 

PCS  3* 

PCS4* 

Target : 


Adj.  R 

.902 

.735 

.908 

.710 

.574 

.403 

MSE 

.007 

.005 

.081 

.256 

.371 

.522 

F** 

122.3 

37.8 

130.9 

33.5 

18.7 

9.9 

KEYTG  1 

.405 

.163 

1.407 

-.107 

-3. 112 

-1.291 

KEYTG  2 

.002 

.046 

-.024 

.539 

-.183 

.544 

KEYTG  3 

-.056 

-.044 

-.153 

.241 

-.661 

.659 

KEYTG  4 

.200 

.057 

.686 

-.506 

1.549 

-1.134 

KEYTG  S 

.351 

.040 

1.445 

2.252 

.606 

-2.332 

KEYTG  6 

.138 

.  188 

.452 

2.409 

1.145 

2.925 

KEYTG  7 

.  104 

.041 

.416 

-1.269 

.688 

.768 

KEYTG  8 

.083 

.048 

.280 

-1.052 

.229 

.551 

Constant 

.026 

.059 

-.799 

-.090 

-.123 

-.073 

Control : 

Adj.  R2 

.646 

.359 

.674 

.214 

.174 

.126 

MSE 

.013 

.004 

.297 

.753 

.822 

.975 

P** 

49.4 

15.8 

55.9 

8.2 

6.6 

4.85 

KEYCN  1 

.123 

.037 

.565 

-.291 

-.044 

.472 

KEYCN  2 

.146 

.021 

.781 

-.347 

.924 

.649 

KEYCN  3 

.258 

.084 

1.443 

-.336 

-1.173 

-.984 

KEYCN  4 

.129 

.059 

.576 

1.689 

-.608 

.841 

Constant 

.127 

.076 

-.588 

-.236 

.110 

.020 

*PCSi  is  the  i^*1  Principal  Component  Score. 


**A11  regressions  were  run  on  all  107  bands  so  the  F  statistic 
is  based  on  8  and  98  degrees  of  freedom  for  the  Target  regres¬ 
sions  and  on  4  and  102  degrees  of  freedom  for  the  Control 
regressions . 
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subsequent  work  that  attempts  to  use  selected  stations  for  the 
summarization  of  precipitation.  The  results  on  the  Control  area 
stations  are  poor  enough  to  indicate  that  the  key  station  ap¬ 
proach  is  not  useful  for  smaller  areas  and  smaller  numbers  of 
stations. 

In  summary,  the  key-station  approach  as  used  here  had  lim¬ 
ited  utility  for  reducing  costs  in  weather  modification  experi¬ 
ments.  However,  for  large  areas,  if  the  mean  is  to  be  used  as 
the  summary  response  measure,  the  approach  could  be  useful,  and 
probably  could  be  improved  on  by  inclusion  of  stations  with  high 
first-component  correlations. 

3.  MULTIVARIATE  ANALYSES  OF 
CLOUD-SEEDING  EFFECTS 

The  models  that  are  covered  below  do  not  exhaust  the  multi¬ 
variate  application  in  the  analysis  of  data  from  cloud-seeding 
experiments.  The  applications  reported  here  ignore  approaches 
that  use  concomitant  variables  because  there  are  indications 
that  the  covariates  from  the  SBA-I  data  were  contaminated  by  the 
seeding  operations  (Bradley,  Srivastava,  and  Lanzdorf;  1977).  In 
a  forthcoming  technical  report  I  have  developed  alternative  meth¬ 
ods  of  analysis  (canonical  correlation  and  multivariate  analysis 
of  covariance)  that  incorporate  the  covariates  to  reduce  error 
variance  in  tests  of  the  effects  of  seeding. 

This  section  considers  two  methods  of  analysis  that  rely  on 
PCA  and  in  each  case,  the  methods  are  applied  to  the  SBA-I  data. 
This  precipitation  data  obviously  violates  the  assumptions  that 
are  necessary  for  parametric  tests  of  significance,  but  some  of 
the  familiar  statistics  are  produced  as  descriptive  measures. 

The  multivariate  precipitation  measurements  could  respond 
to  seeding  in  several  ways.  Cloud-seeding  could  have  no  effect, 
it  could  affect  the  means  of  the  vectors  of  precipitation  measure¬ 
ments  while  leaving  the  covariances  between  stations  unchanged, 
it  could  leave  the  means  unaffected  but  have  an  effect  on  the 
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covariance  structure,  or  both  means  and  covariances  could  be 
affected  by  seeding.  If  seeding  has  no  effect,  then  the  princi¬ 
pal  components,  and  vectors  of  principal  component  scores,  esti¬ 
mated  from  all  available  observations  should  be  essentially  the 
same  for  the  seeded  and  unseeded  experimental  units. 

If  however,  seeding  affects  the  means  and/or  covariances, 
the  analysis  should  consider  summary  measures  that  axe  estimated 
separately  for  the  seeded  and  nonseeded  observations.  The  two 
subsections  below  present  analyses  based  on  these  alternative 
assumptions. 

4.1  Analysis  of  Principal  Components  from  Pooled-Data  Covariance 
Matrices 

The  discussion  and  results  in  Section  3  implicitly  assume 
that  the  means  and  the  covariance  structures  are  identical  for 
seeded  observations  and  unseeded  observations.  If  this  is  true, 
the  means  of  the  principal  component  score  vectors  should  be  the 
same  for  seeded  and  unseeded  observations.  Thus  a  simple  test 
for  the  effects  of  seeding  would  be  to  test  the  hypothesis  that 
the  mean  vector  of  principal  component  scores  for  seeded  obser¬ 
vations  equals  the  mean  vector  of  principal  component  scores  for 
unseeded  observations. 

It  is  possible  that  scores  for  some  components  have  differ¬ 
ent  means  although  the  vectors  of  means  do  not  differ  significant¬ 
ly.  For  example,  univariate  tests  of  the  orographic  component 
scores  could  indicate  differences  for  seeded  and  unseeded  bands 
even  though  a  vector  test  could  obscure  a  significant  difference. 
The  subsections  below  describe  methods  for  the  analysis  of  vec¬ 
tors  of  scores  and  for  the  analysis  of  univariate  scores  as  well. 


4.1.1  A  Multivariate  Comparison  of  Principal  Component  Scores 
2 

Hotelling's  T  statistic  is  an  appropriate  test  statistic 
for  testing  the  hypothesis  that  jjj  =  ^2  when  the  two  populations 
of  interest  have  multivariate  normal  distributions  with  a  common 


covariance  matrix,  £,  of  full  rank  (the  location  vectors  may  dif¬ 
fer,  however).  These  assumptions  do  not  hold  for  the  principal 
components  scores  of  precipitation  data,  especially  when  the 
principal  components  are  extracted  from  the  correlation  matrix 
(see  Morrison  (1976)  for  a  discussion  of  the  distributional  prob¬ 
lems  here).  Nonetheless,  if  these  limitations  are  recognized  as 

2 

one  interprets  the  results.  Hotelling's  T  statistic  can  indicate 
similarities  in  the  vectors  of  mean  principal  component  scores, 
since  a  true  null  hypothesis  would  not  be  rejected  although  the 
power  to  reject  a  false  null  hypothesis  is  limited. 

For  the  principal  component  results  given  in  Section  3, 

2 

Hotelling's  T  was  computed  to  compare  the  vectors  of  the  means 

of  the  first  three  principal  component  scores  for  the  Target 

2 

Area,  the  computed  T  value  corresponded  to  an  F-ratio  of  1.016 
with  3  and  103  degTee  of  freedom  (At  a  =  .05  the  critical  value 
of  F  is  2.7  with  3  and  100  degrees  of  freedom).  If  the  distri¬ 
butional  assumptions  held,  this  result  would  indicate  negligible 
effects  for  seeding. 

4.1.2  Univariate  Comparisons  of  Principal  Component  Scores 

As  the  introduction  to  this  section  notes,  univariate  com¬ 
parisons  could  indicate  differences  for  the  means  of  the  indivi¬ 
dual  components  that  get  "washed  out"  in  the  multivariate  statis¬ 
tics.  Of  course  if  a  series  of  univariate  tests  of  significance 
are  produced,  such  multiple-comparison  tests  suffer  from  a  loss 
of  power.  In  the  present  case,  the  data  do  not  meet  the  usual 
assumptions  for  tests  of  significance,  so  the  statistics  are 
produced  for  descriptive  purposes  only. 

The  t-statistics  (the  difference  in  the  estimated  mean 
principal  component  scores  for  seeded  and  nonseeded  data,  divi¬ 
ded  by  the  pooled  standard  deviations  of  the  principal  component 
scores)  for  the  principal  component  scores  of  the  Target  Area 
are  given  in  Table  5.  Only  the  first  component  means  are  very 
different  and  this  indicates  that  if  there  are  differences  in 
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the  principal  components  from  the  pooled  data,  the  differences 
seem  to  be  associated  with  the  band  means. 


TABLE  -5:  UNIVARIATE  T-STATISTICS  FOR  THE  FIRST  THREE 
PRINCIPAL  COMPONENT  SCORES  FROM  THE 
POOLED  TARGET  AREA  DATA 


Component 


t-Statistics  (N  =  107  bands) 


1 

2 

3 


1.67 

.52 


.07 


This  concludes  comparisons  of  the  seeded  and  nonseeded  sub¬ 
sets  of  components  based  on  a  pooled-data  set.  The  next  subsec¬ 
tion  presents  results  obtained  when  the  data  was  split  into 
seeded  and  unseeded  bands  to  produce  two  sets  of  principal  com¬ 
ponent  estimates. 

4.2  The  Effects  of  Seeding  on  the  Covariance  Structure 

Separate  PCA  estimates  (one  from  the  seeded  experimental 
units  and  one  from  the  unseeded  experimental  units)  can  provide 
insight  into  the  effects,  if  any,  of  cloud-seeding  on  the  covari¬ 
ance  structure  of  the  precipitation  measurements.  First,  the 
data  from  the  stations  are  separated,  then  standardized  in  the 
manner  of  Section  (2.1)  to  produce  separate  estimates  of  correl¬ 
ation  matrices  for  the  seeded  and  unseeded  observations.  The 
principal  components  estimated  from  these  correlation  matrices 
should  reflect  changes,  if  any,  in  the  station-to-station  covari¬ 
ances  that  are  produced  by  seeding.  Tests  for  significant  dif¬ 
ferences  are  again  limited  since  the  data  does  not  meet  the  req¬ 
uisite  assumptions  for  tests  of  significant  differences. 

A  summary  of  the  results  on  the  estimated  eigenvalues  for 
the  seeded  and  unseeded  data  sets  as  well  as  the  previously- 
reported  results  on  the  combined  data  set  are  given  in  Table  6. 
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This  table  shows  that  the  first  seeded  vs.  nonseeded  eigenvalues 
differ  by  6%  of  the  total  variance  for  the  Target  Area  versus 
1.6%  the  total  variance  for  the  Control  Area.  Even  when  the 
requisite  assumption  of  multivariate  normality  holds,  the  vari¬ 
ances  of  the  sampling  distributions  are  so  large  that  the  Target 
Area  difference  would  be  a  nonsignificant  difference  (Press  1972). 
The  results  do  indicate  that  the  seeding  operations  produced  pre¬ 
cipitation  measurements  over  the  Target  Area  raingage  network 
that  were  more  correlated  with  the  first  component  (and  thus, 
the  band  means)  of  the  precipitation  and  further,  the  measurements 
were  less  correlated  with  the  second  (orographic)  component.  The 
same  result  does  not  hold  for  the  control  area  data  where  the 
differences  were  less  marked  for  all  components. 

The  changes  in  the  first  two  eigenvalues  could  be  accompanied 
by  changes  in  the  patterns  of  station-component  correlations  (as 
in  Figure  1-6).  Although  graphs  are  not  reproduced  here  because 
of  space  considerations,  the  basic  patterns  of  signs  for  the 
first  three  components  were  similar  to  the  patterns  of  signs  in 
Figures  1-6  for  both  the  seeded  and  nonseeded  component  estimates. 
Of  course,  the  relative  heights  of  the  spikes  (which  represent 
the  stations-to-component  correlations)  changed  to  reflect  in¬ 
creases  or  decreases  in  correlations. 

In  conclusion,  the  proposed  analysis  and  the  results  obtained 
with  SBA-I  data  suggest  that  analyses  of  weather  modification 
data  should  consider  the  effects  of  seeding  on  the  covariance 
structure  of  the  precipitation  measurements.  The  data  here  rep¬ 
resent  such  gross  violations  of  the  requisite  assumptions  for 
inference  that  no  conclusions  can  be  reached  on  the  effects  of 
seeding  on  the  covariance  structure  of  precipitation  measurements. 
In  general,  the  nature  of  the  distributions  of  precipitation 
measurements  precludes  tests  based  on  assumptions  of  normality. 
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5.  CONCLUSIONS  AND  IMPLICATIONS 


As  is  the  case  with  many  other  analyses  of  weather  modifica¬ 
tion  data,  the  conclusions  are  limited  and  the  implications  for 
further  analysis  nearly  boundless.  The  premise  that  multivariate 
models  can  be  applied  to  data  from  weather  modification  data  is 
supported  by  the  results  here. 

The  applications  of  PCA  to  precipitation  measurements  produced 
readily- interpretable  principal  components.  For  both  Control  and 
Target  Areas,  the  first  components  were  almost  perfectly  correlated 
with  the  simple  average  of  the  experimental  units,  the  second  component 
was  readily  interpreted  as  an  orographic  component,  while  the  third 
component  was  identified  as  an  east-west  component.  The  use  of  PCA 
to  select  a  subset  of  key  stations  has  utility  for  large  areas  but 
the  method  produced  inadequate  summary  measures  for  small  areas  with 
few  stations. 

The  use  of  multivariate  methods  to  analyze  the  effects  of  cloud- 
seeding  is  also  limited.  The  data,  particularly  the  standardized  co- 
variance  matrix  (i.e.  the  correlation  matrix)  fails  to  meet  the  multivariate- 
normality  assumption  so  that  the  results  do  not  allow  the  usual  inferences 
as  to  the  effects  of  cloud-seeding.  The  estimates  obtained  from  the  SBA-I 
data  indicate  that  the  seeded  bands  produced  precipitation  measurements 
that  are  uniformly  more  highly  correlated  with  the  band  means  than  the 
measurements  for  the  unseeded  bands. 
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APPENDIX 


This  Appendix  provides  detailed  information  from  the 
principal-components  results.  The  results  are  organized  by 
Target  Area  and  Control  Area  as  well  as  whether  the  estimates 
were  based  on  all  bands,  nonseeded  bands  only,  or  seeded  bands 
only.  The  first  Tables  give  descriptions  of  the  stations  and 
the  other  tables  give:  means  and  standard  deviations;  communal- 
ities  and  eigenvalues  for  the  principal  components;  station-to- 
component  correlations;  as  well  as  the  station- to- station  correla¬ 
tion  matrix  (For  the  correlation  matrices,  in  the  interest  of 
space,  only  the  correlations  that  are  based  on  all  bands  are 
given.  Copies  of  the  nonseeded-and  seeded-band  correlation 
matrices  are  available  on  request) . 
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Table  A.l.T.  TARGET  AREA  STATIONS  USED 

IN  THE  PRINCIPAL  COMPONENT  ANALYSIS 


STATION 

NAME 

NO.  MISSING 
BANDS 

ORIG.  STAT. 
NUMBER 

LATITUDE 

(DEG..01S) 

LONGITUDE 

(DEG..01S) 

ALTITUDE 

(METERS) 

E1253 

7 

25 

34.58 

119.98 

238 

E1682 

0 

26 

34.23 

118.62 

276 

E3751 

0 

29 

34.27 

118.40 

335 

E5756 

20 

36 

35.05 

118.17 

’.834 

E7735 

8 

40 

34.75 

118.73 

1377 

E8752 

4 

47 

35.15 

119.47 

312 

E8832 

20 

48 

35.13 

118.43 

1207 

LOOIO 

22 

55 

34.08 

118.45 

178 

LOOll 

6 

56 

34.12 

118.42 

264 

L0033 

1 

57 

34.33 

118.40 

457 

L0046 

4 

58 

34.30 

118.18 

706 

L0053 

0 

59 

34.30 

118.12 

1104 

L0054 

22 

60 

34.35 

118.05 

1319 

L0128 

18 

62 

34.60 

118.57 

633 

L0179 

4 

64 

34.17 

118.07 

360 

L0191 

4 

65 

34.07 

118.20 

122 

L0213 

1 

68 

34.07 

118.35 

53 

L0250 

15 

69 

34.45 

118.20 

8C0 

L0259 

2 

70 

34.28 

118.60 

389 

L0292 

13 

72 

34.15 

118.52 

328 

L0303 

1 

73 

34.13 

118.12 

244 

L0372 

3 

75 

34.53 

118.52 

482 

L0373 

7 

76 

34.23 

118.22 

678 

L0434 

13 

77 

34.13 

118.75 

244 

L0435 

10 

78 

34.08 

118.70 

183 

L0438 

6 

79 

34.02 

118.80 

35 

L0466 

5 

80 

34.35 

118.35 

982 

L0694 

5 

83 

34.28 

118.28 

465 

L1014 

3 

86 

34.00 

118.10 

52 

L1017 

12 

87 

34.48 

118.02 

1000 

L1074 

17 

88 

34.38 

118.15 

1707 

L1107 

9 

89 

34.23 

118.33 

354 

L1181 

18 

91 

34.33 

118.45 

1495 

H  124 

4 

98 

34.58 

118.37 

930 

M.  223 

12 

99 

34.42 

119.68 

2 

M  226 

10 

100 

34.45 

119.95 

34 

11  228 

1 

101 

34.45 

119.68 

262 

M  230 

2 

102 

34.52 

119.68 

473 

N  13 

4 

112 

34.53 

119.95 

1220 

N  14 

4 

113 

34.63 

119.87 

348 

R  13 

9 

114 

34.55 

119.87 

282 

N  A 

19 

120 

34.53 

120.02 

1067 
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Table  A.l.T.  (Continued) 


STATION 

NAME 

NO.  HISSING 
BANDS 

ORIG.  STAT. 
NUMBER 

LATITUDE 

(DEG..01S) 

LONGITUDE 
(DEG,. 018) 

ALTITUDE 

(METERS) 

H  B 

3 

121 

34.43 

119.83 

3 

S  208 

10 

128 

34.40 

119.52 

3 

S  210 

4 

129 

34.45 

119.65 

168 

S  211 

3 

130 

34.45 

119.78 

122 

S  221 

3 

134 

34.98 

119.67 

662 

S  231 

9 

135 

34.45 

119.57 

427 

S  232 

2 

136 

34.48 

119.50 

633 

S  234 

3 

138 

34.42 

119.70 

24 

S  238 

10 

141 

34.78 

119.65 

1524 

S  242 

19 

142 

64.48 

119.80 

305 

V  085 

21 

148 

34.38 

119.23 

335 

V  138 

14 

150 

34.63 

119.43 

1479 

V  165 

4 

151 

34.47 

119.25 

457 

V  166 

3 

152 

34.42 

119.35 

229 

V  167 

0 

153 

34.28 

119.27 

91 

V  168 

2 

154 

34.20 

118.20 

11 

V  171 

9 

157 

34.40 

118.88 

137 

V  172 

13 

158 

34.52 

118.77 

351 

V  173 

0 

159 

34.43 

119.08 

305 

V  174 

20 

160 

34.68 

119.35 

1091 

V  190 

4 

164 

34.28 

119.00 

168 

V  191 

0 

165 

34.32 

118.88 

183 

V  192 

22 

166 

34.25 

118.85 

168 

V  194 

6 

168 

34.20 

119.02 

37 

V  198 

16 

170 

34.57 

119.17 

1055 

V  206 

4 

171 

34.32 

118.97 

183 

V  207 

18 

172 

34.50 

119.38 

500 

V  209 

6 

173 

34.73 

119.10 

1570 

V  221 

1 

174 

34.35 

119.42 

3 

V  225 

4 

176 

34.38 

119.15 

213 

Table  A.l.C.  CONTROL  AREA  STATIONS  USED 

IN  THE  PRINCIPAL  COMPONENT  ANALYSIS 


STATION 

NAME 

NO.  MISSING 
BANDS 

ORIG.  STAT. 
NUMBER 

LATITUDE 
(DEG, .01S) 

LONGITUDE 

(DEG..01S) 

ALTITUDE 

(METERS) 

E4144 

10 

30 

35.10 

120.38 

218 

E7946 

0 

43 

34.90 

120.45 

72 

N  01 

12 

103 

34.83 

120.53 

91 

N  02 

10 

104 

34.62 

120.47 

98 

N  03 

0 

105 

34.75 

120.28 

177 

N  04 

4 

106 

34.60 

120.20 

104 

N  05 

5 

107 

34.65 

120.12 

238 

N  06 

19 

108 

34.73 

120.08 

375 

S  201 

0 

122 

34.78 

120.33 

195 

S  202 

8 

123 

34.73 

120.23 

220 

S  203 

22 

124 

34.73 

120.37 

305 

S  204 

5 

125 

34.75 

120.28 

172 

S  205 

18 

126 

34.68 

120.45 

91 

S  206 

7 

127 

34.47 

120.23 

9 

S  217 

18 

133 

34.50 

120.48 

6 

S  233 

0 

137 

34.62 

120.20 

114 

S  235 

0 

139 

34.95 

120.45 

61 

S  236 

1 

140 

34.53 

120.18 

201 

S  251 

13 

144 

34.58 

120.50 

341 

Table  A.2.TA.  Means  and  Standard  Deviations  for  All  Bands 
Target  Area. 


VARIABLE 

mean 

STANDARD  DEV 

CASES 

T90I 

.2815 

.3*16 

100 

T  002 

.1838 

.2368 

107 

T003 

.203* 

.2987 

107 

TOO** 

.05?! 

.1180 

87 

T005 

.1316 

.1982 

99 

T006 

.05g2 

.09*0 

103 

TOOT 

.082* 

.1185 

07 

toob 

.1TT9 

•  26T5 

85 

T009 

.2135 

.3292 

101 

T0I0 

.2*39 

.3*11 

106 

T0I1 

.3131 

.*890 

103 

T012 

.2926 

.*716 

107 

T013 

.1958 

.28*1 

85 

T01* 

.2521 

.3399 

89 

TOJS 

.2388 

.3902 

103 

T016 

.!Se3 

.2175 

103 

TOIT 

.1728 

.25*8 

106 

T01B 

.1225 

.1987 

92 

T019 

.1952 

.260* 

105 

T02O 

.2228 

.3819 

9* 

T02I 

.2153 

.3039 

106 

T022 

.1795 

•  ZS*3 

10* 

T023 

.3025 

.*283 

100 

T02*> 

.2055 

.3753 

9* 

T?25 

.2*97 

.3936 

97 

T026 

.1972 

.2870 

101 

T02T 

.2366 

.3135 

102 

T028 

.2136 

.3251 

102 

T029 

.15i5 

.1986 

10* 

TO  30 

.1002 

.  1558 

95 

T03I 

.253* 

.3308 

90 

TO  32 

.2026 

.3328 

98 

T033 

.2176 

.2831 

89 

TO  3* 

.178* 

.2776 

103 

TO  35 

.25o2 

.2988 

95 

1  V 

T03T 

" 

.  .2S6T 
.3*11 

9/ 

106 

T038 

.02*9 

.5730 

105 

T039 

.*0|8 

.*5*5 

103 

T  0*0 

.25?* 

.2385 

103 

TO*  I 

.3642 

.*292 

98 

T0*2 

.3606 

.*167 

88 

T0«3 

.2306 

.26*3 

10* 

TO** 

.2*01 

.2996 

97 

T0*5 

.3096 

.38S* 

103 

T0*6 

•  25|  0 

.2962 

10* 

T0*T 

.0689 

.0929 

10* 

TO*B 

«3S|8 

.*189 

98 

T  0*9 

.5168 

.7517 

105 

T050 

.2*13 

.3017 

10* 

T051 

.I838 

.2522 

97 

T  052 

.33*9 

.3580 

88 

T053 

.323* 

.3901 

86 

T05* 

T055 

‘Alii 

.2*58 

.3613 

93 

103 

T056 

.2999 

.3558 

10* 

T05T 

.203* 

.2793 

107 

T058 

.1903 

.2767 

105 

T059 

.2378 

.2825 

98 

T060 

.19*0 

.2352 

9* 

T061 

.3086 

.*1*8 

107 

T  062 

.1859 

.2638 

87 

T063 

.208* 

.2892 

103 

T06* 

.2052 

.2S36 

107 

T065 

.1796 

.2**3 

85 

T066 

.1672 

.2*95 

101 

TOOT 

.308* 

.*0*2 

91 

TOOB 

.2038 

.255* 

103 

T009 

.3631 

.*387 

89 

TOTO 

.135S 

.1831 

101 

TOTl 

.1991 

.2*67 

106 

T072 

.3003 

.3911 

103 

36 


Table  A.2.TN 


Means  and  Standard  Deviations  for  Nonseeded  Bands 
Target  Area. 


VARIABLE 

HCaN 

STANDARD  dev 

CASES 

T001 

.2270 

.2858 

47 

T002 

•  IL?5 

.1884 

51 

TOO) 

•  1490 

.2179 

51 

TOOL 

.0576 

.1358 

41 

T005 

.1158 

.1910 

48 

T006 

.0L?9 

.0922 

49 

TOOT 

.0633 

.0915 

39 

T008 

.1183 

.1864 

41 

T009 

.1525 

.2414 

48 

T010 

.1988 

.30S1 

51 

T0I1 

•  2L0O 

.4171 

50 

TO  12 

.2106 

.3767 

51 

TOD 

.1266 

.2015 

38 

TOIL 

•  19|L 

.2844 

42 

T01S 

.2008 

.3112 

50 

T016 

.1208 

.1633 

48 

T0I7 

.1451 

.2112 

51 

T018 

.103L 

.1790 

44 

TOD 

.1368 

.1807 

50 

T  020 

.1947 

.3690 

45 

T021 

.1737 

.2437 

51 

T022 

.1296 

.1818 

50 

T023 

.2373 

.3563 

48 

T02L 

.2360 

.3510 

45 

TC25 

.1983 

.3022 

48 

T026 

.1913 

.2914 

48 

T027 

.1878 

.2708 

49 

T028 

.1835 

.2745 

49 

TC29 

.1229 

.1596 

51 

T030 

.0798 

.U43 

45 

T031 

.1974 

.2716 

43 

T032 

.1638 

.2988 

47 

T033 

.1550 

.1759 

44 

T03L 

.1176 

.1737 

50 

T035 

.2187 

.2446 

46 

T036 

.1946 

.1862 

48 

T037 

.2557 

.2903 

51 

T038 

.3204 

.4140 

SO 

T039 

.3244 

.3622 

48 

TOLO 

.2037 

.2110 

49 

TOLI 

.3036 

.3899 

45 

T0L2 

.2760 

.3087 

40 

TOO 

.1824 

.1857 

50 

TOLL 

.2311 

.3142 

46 

TOLS 

.2490 

.3101 

49 

T0L6 

.2037 

.2208 

49 

T0L7 

.0553 

.0856 

49 

TOLS 

.2860 

.3347 

47 

T0L9 

.4120 

.6411 

50 

TOSO 

.2028 

.2572 

50 

T051 

.1467 

.  16S2 

45 

T  052 

.2756 

.2876 

45 

T0S3 

.2931 

.3952 

3~ 

TOSL 

.1591 

.1857 

43 

TOSS 

.2610 

•  3524 

48 

TCS6 

.2708 

.3573 

50 

T0S7 

.1802 

.2649 

51 

T05L 

.1591 

.1857 

43 

T0S5 

.2610 

.3524 

48 

T056 

.2708 

.3573 

50 

T057 

.1802 

.2649 

51 

TOSS 

.1726 

.2942 

50 

T0S9 

.1915 

•  2S62 

48 

T060 

.1307 

.1920 

45 

TC61 

.2455 

.3577 

51 

T062 

.1578 

.2070 

40 

7063 

.1762 

.2626 

48 

T06L 

.1641 

.2144 

51 

T  065 

.1707 

.2357 

42 

T066 

.1479 

.2637 

48 

TC67 

.2267 

.2954 

43 

T068 

.1508 

.1969 

49 

T  069 

.2909 

.3829 

43 

T070 

.1172 

.1873 

46 

TOTI 

.1720 

.1973 

51 

T072 

.2386 

.3470 

49 

Table  A 


2.TS.  Means  and  Standard  Deviations  for  Seeded  Bands 
Target  Area. 


VARIABLE 

MEAN 

STANOARO  OEV 

CASES 

T001 

.3208 

.3806 

53 

T002 

.221* 

.2698 

56 

TOOJ 

.2529 

.3516 

56 

TOO* 

.0567 

.1011 

*6 

TOOS 

.1*65 

.2055 

51 

T006 

.0720 

.09** 

5* 

TOOT 

.0979 

.1355 

*8 

TOOS 

.233* 

.3176 

** 

T  009 

•  26«7 

.3863 

53 

T010 

.28$6 

.3693 

55 

1011 

.3705 

.5*51 

S3 

T0I2 

.367  3 

.536* 

56 

T013 

.2S|7 

.3279 

*7 

T01* 

.306* 

.3776 

*7 

T015 

.27*7 

.*525 

S3 

T016 

.1909 

.2526 

55 

TOIT 

.»9fl5 

.2891 

55 

T01S 

.1*00 

.2155 

*8 

TO  IS 

.2*8* 

.3081 

55 

T020 

.2*86 

.3955 

*9 

T021 

.2S38 

.3*85 

55 

T  022 

.2257 

.3010 

5* 

T023 

.3627 

.*810 

52 

102* 

.2927 

.3980 

*9 

1025 

.3000 

.*638 

*9 

T026 

.2026 

.2856 

53 

T02T 

.2817 

.3**7 

53 

T  028 

.2*15 

.3662 

S3 

T029 

.1791 

.2282 

53 

T030 

.1186 

.1722 

SO 

T03I 

.30*7 

.3725 

*7 

TO  32 

.2382 

.3606 

51 

TO  33 

.2789 

.3*98 

*5 

T03* 

.2358 

.3*0* 

53 

T035 

.2971 

.3*02 

*9 

TO  36 

.2761 

.3073 

*9 

T037 

.3**5 

.3800 

55 

T  038 

.5198 

.676* 

55 

TO  39 

.9695 

•  S1S8 

55 

TO  AO 

.2967 

.25*8 

5* 

TOM 

.*157 

.*57* 

53 

T  0*2 

.9*58 

.*785 

48 

T0*3 

•  28*9 

.3139 

5* 

TO** 

.2*82 

.2887 

51 

T0*5 

.  36*6 

.*38* 

5* 

T  0*6 

.2931 

.3*66 

55 

T0*7 

'  .0811 

(  .0981 

55 

TO*8 

.9125 

.*791 

SI 

T0*\ 

.6120 

.83*2 

55 

T  050 

.2770 

.3362 

5* 

T051 

.2160 

.306* 

52 

052 

.3970 

.*138  . 

*3 

T053 

.3*85 

.3882 

*7 

TOS* 

.2206 

.2865 

50 

T055 

.3195 

.3700 

55 

T056 

.3269 

.3556 

5* 

T057 

.22*5 

.2927 

56 

TOSS 

.2178 

.2605 

55 

T059 

.2822 

.301* 

50 

T060 

.2522 

.2571 

*9 

T06I 

.3661 

.*563 

56 

T062 

.2098 

.30*1 

*7 

T063 

.2365 

.3102 

55 

T06* 

.2*27 

.281* 

56 

1065 

.18** 

.25*9 

*3 

T066 

.18*7 

.2371 

53 

T067 

.38 15 

.*726 

*8 

T068 

.2S|9 

.292* 

5* 

T  069 

.9307 

.*795 

*6 

T070 

.1509 

.1797 

55 

T0T1 

.22*2 

.28*6 

55 

T072 

.3563 

.*225 

5* 
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TABLE  A.2.C.  Means  and  Standard  Deviations  for  All  Bands,  Non- seeded 
Bands,  and  Seeded  Bands:  Control  Area. 


Table  A.2.C 


Means  and  Standard  Deviations  for  All  Bands,  Non 
seeded  Bands,  and  Seeded  Bands:  Control  Area. 
(Continued) . 


NOT  SEEDED  BANDS 


VARIABLE 

rean 

stanoaro  oev 

CASES 

C01 

.223* 

.2571 

*7 

C02 

.1753 

.207* 

51 

CO  3 

.1630 

.1761 

*6 

CCA 

•  222* 

.1876 

*5 

COS 

.10*0 

.1-731 

51 

C86.. 

-  -»?.0?2_ 

.  .2002 

.51 

COT 

•  ft  122 

.1669 

*v 

COO 

.2086 

.2030 

A* 

coo 

.1971 

.1617 

51 

CIO 

•  207A 

.1686 

*7 

Cll 

.2A5S 

.1018 

*0 

CI2 

.2208 

.1033 

*0 

C13 

.<80S 

.1**2 

*1 

CIA 

.10*0 

.1705 

*7 

CIS 

.23*5 

.1087 

A* 

C16 

.22|0 

.2170 

SI 

Cl  7 

.1790 

.1850 

51 

CIS 

.28*0 

.3007 

50 

CIO 

.2S|  | 

.20*0 

*5 

SEEDED  BANDS 


variable 

HtAN 

STANOARO  OEV 

CASES 

COI 

.2620 

.2380 

SO 

C02 

,1737 

.1017 

56 

COS 

'  .  .1890 

.2262 

69 

CO* 

.2787 

.2865 

52 

COS 

.2075 

.1013 

56 

C06 

.2167 

.23*0 

52 

COT 

.1057 

.1861 

53 

COB 

.2380 

.220* 

64 

COO 

.2121 

.1009 

56 

CIO 

.2320 

.2108 

52 

Cll 

.2*13 

.1808 

45 

C12 

.2275 

.2099 

53 

CIS  . 

.10]S 

.1676 

48 

Cl* 

.2521 

.2608 

S3 

CIS 

.32*0 

.3*26 

45 

C16 

.2*05 

.2*62 

56 

Cl  7 

.178* 

.18*8 

56 

■?C18 

.0112 

.*612 

56 

CIO 

.3251 

.3005 

49 

40 


Component 

EIGENVALUE 

pci  or  »*« 

CUM  PCT 

i 

51*33061 

71.3 

71.3 

2 

4*60169 

6.7 

78.0 

3 

4.23604 

5.9 

83.8 

4 

2.07392 

2.9 

86.7 

S 

1* 49353 

1.9 

88.7 

6 

.99920 

1.4 

90.0 

7 

iioa'ro 

1.1 

91.2 

a 

.13992 

1.0 

92.2 

9 

.58893 

.8 

93.0 

19 

.55613 

.8 

93.8 

11 

.54899 

.8 

94.6 

12 

.51317 

.7 

95.3 

13 

.45599 

.6 

95.9 

14 

.408/3 

.6 

96.5 

IS 

.40S98 

.6 

97.0 

16 

.31322 

.4 

97.5 

1? 

.30512 

.4 

97.9 

IS 

.27041 

.4 

98.3 

19 

•247i2 

.3 

98.6 

20 

.21995 

.3 

96.9 

21 

.21186 

.3 

99.2 

22 

.18759 

.3 

99.5 

23 

.1579/ 

.2 

99.7 

24 

.15614 

.2 

99.9 

25 

.14442 

.2 

100.1 

26 

.13096 

.2 

100.3 

27 

.12569 

.2 

100.5 

28 

.11424 

.2 

100.6 

29 

.10943 

.2 

100. a 

30 

•09550 

.1 

100.9 

31 

.08753 

.1 

101.0 

32 

.08524 

.1 

101.1 

33 

.07719 

.1 

101.3 

34 

.07143 

.1 

101.4 

35  ... 

.06195 

.  .1 

101.4 

36 

•OS601 

.1 

101.5 

37 

.05203 

.1 

101.4 

38 

.04521 

.1 

101.7 

39 

.03861 

.1 

101.7 

40 

.03418 

.0 

101.8 

41 

.031*2 

.0 

101.8 

42 

.02751 

.0 

101.8 

43 

.02*53 

.0 

101.9 

44 

.02165 

.0 

101.9 

45 

.02046 

.0 

101.9 

46 

.01782 

.0 

102.0 

47 

.01548 

.0 

102.0 

48 

.01181 

.0 

102.0 

49 

.01005 

.0 

102.0 

SO 

.00844 

.0 

i02.0 

51 

.00552 

.0 

102.0 

52 

.00338 

.0 

102.0 

53 

.00298 

.0 

102.0 

54 

••00048 

-.0 

102.0 

55 

-.00330 

-.0 

102.0 

56 

-.00858 

-.0 

102.0 

57 

-.01065 

-.0 

102.0 

58 

-.91356 

-.0 

102.0 

59 

-.01967 

-.0 

102.0 

60 

-.02369 

-.0 

101.9 

61 

-.02946 

-.0 

101.9 

62 

-.03364 

-.0 

101.8 

63 

-.03858 

-.1 

101.8 

64 

-.042/1 

—  i 

101.7 

65 

-.06526 

-.1 

101.8 

66 

-.07880 

-.1 

181.5 

67 

-.09325 

-.1 

101.4 

68 

-.12161 

-.2 

101.2 

69 

-.16650 

-.2 

101.0 

70 

— .i9940 

-.3 

100.7 

71 

-.25129 

-.3 

100.4 

72 

-.26664 

-.4 

100.0 
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Table  A.3.TN.  Communalities  and  Eigenvalues  for  Nonseeded  Bands: 
Target  Area. 


VARIABLE 

EST  COMMUNAL! TV 

Component 

EIGENVALUE 

PCI  or  VAR 

CUH  PCT 

root 

.939*0 

i 

*0. 7 I9VV 

67.7 

67.7 

1002 

.9690* 

2 

7.36707 

10.2 

77.9 

T003 

.90888 

3 

*.31830 

6.0 

83.9 

TOO* 

.91913 

* 

2.65997 

3.7 

87.6 

1005 

•07?IT 

5 

1  •  vsiba 

2.7 

90.3 

T  006 

.89271 

6 

1.6*288 

2.1 

92.6 

T007 

.01226 

7 

1.19896 

1.7 

9*. 2 

T008 

.958*7 

8 

.V1036 

1.3 

95.5 

7009 

.958*7 

9 

.86993 

1.2 

96.7 

T0I0 

.0012* 

10 

.70233 

1.9 

77.7 

TOU 

.936** 

11 

.608*6 

.8 

98.5 

7012 

.95163 

12 

.59521 

.8 

99.* 

T013 

.07*7* 

13 

.53239 

.7 

100.1 

TOI* 

.9*02“ 

1* 

.*5697 

.6 

ioo.7 

T0IS 

.9**77 

IS 

.37920 

.5 

101.3 

T016 

.9**71 

16 

.33205 

.5 

101.7 

TOI  T 

.9**08 

17 

.30013 

.* 

102.1 

T018 

.9*6S5 

18 

.2S201 

.* 

102.5 

7019 

.9680* 

19 

.22079 

.3 

102.0 

T  020 

.96771 

20 

.19988 

.3 

103.1 

7021 

.9**71 

21 

.18*27 

.3 

103.3 

TQ22 

.92508 

22 

.1*562 

.2 

103.5 

T023 

.92697 

23 

.13320 

.2 

103.7 

T02* 

.96771 

2* 

•  10863 

.2 

103.9 

T025 

.95380 

25 

.10083 

.1 

10*. 0 

7026 

.95380 

26 

.09562 

.1 

10*. 1 

T027 

.91765 

27 

.08*88 

.1 

10*.3 

T028 

.96016 

28 

.07311 

.1 

10*.* 

T029 

.08657 

29 

.070l8 

.1 

10*. 5 

T030 

.9*655 

30 

.06373 

.1 

i0*.6 

T031 

.95163 

31 

.06011 

.1 

10*. 6 

TO  32 

.96016 

32 

.0551* 

.1 

10*. 7 

T033 

.90888 

33 

.0*79* 

.1 

10*. 8 

TO  3* 

.92090 

3* 

.03703 

.1 

10*. 8 

T03S 

.9*1*6 

35 

.0321* 

.0 

10*. 9 

.93t*8 

36 

.0300* 

.0 

10*. 9 

T03T 

.9527* 

37 

.02797 

.0 

109.0 

T038 

.939*6 

38 

.02270 

.0 

105.0 

T039 

.97331 

39 

.02030 

.0 

105.0 

T0*0 

.90291 

*0 

.01669 

.0 

105.0 

T0*1 

.93*19 

*1 

.01223 

.0 

10S.1 

T0*2 

.97331 

*2 

.01122 

.0 

ios.i 

T0*3 

.9*607 

*3 

.00927 

.0 

105.1 

TO** 

•86857 

** 

.00799 

.0 

105.1 

T0*5 

.97i3* 

*5 

.00*12 

.0 

Ios.i 

T  0*8 

.95*52 

*6 

.00016 

.0 

105.1 

T0*7 

.09271 

*7 

-.000*7 

-.0 

ios.i 

T0*8 

.92293 

*8 

-.0033* 

-.0 

105.1 

T0*9 

.93227 

*9 

-.00855 

-.0 

105.1 

_ T650 

.  .97.3* 

■  A 

-.00976 

-.0 

105.1 

TOSI 

•09 j *9 

51 

-.01333 

-.0 

105.1 

T052 

.95*52 

52 

-.01*0? 

-.0 

105.0 

T053 

.93678 

S3 

-.01593 

-.0 

ios.o 

T05* 

.9039* 

5* 

-.01962 

-.0 

ios.o 

T055 

.938*9 

55 

-.02687 

-.0 

10*. 9 

T056 

.938*9 

56 

-.03599 

-.0 

10*.9 

T057 

.06*36 

57 

-.03996 

-.1 

to*. 8 

T058 

.89678 

58 

-.0*620 

-.1 

10*. 8 

T059 

.91*32 

59 

-.05060 

-.1 

10*.  7 

T060 

.0*629 

60 

-.07506 

-.1 

10*.6 

T06I 

.9*189 

61 

-.07832 

-.1 

10*. S 

T  062 

.9039* 

62 

-.08927 

-.1 

10*.* 

T063 

.91763 

63 

— • 10920 

-.2 

10*. 2 

TOO* 

.92528 

6* 

-.12236 

-.2 

10*. 0 

T06S 

.9*3*8 

65 

-.162*8 

-.2 

103.8 

T066 

.089*3 

66 

- .20*93 

-.3 

103.5 

T067 

.91051 

67 

-.28138 

103.1 

T068 

.9*3*8 

68 

-.30362 

-.* 

102.7 

T069 

.92782 

69 

- • J*6*8 

-.5 

102.2 

TOTO 

.002*1 

70 

-.*5685 

-.6 

101.6 

T07I 

.06392 

71 

-.52091 

-.7 

100. 9 

T072 

.9*189 

72 

-.637*8 

-.9 

100.0 

42 


Table  A.3.TS.  Communalities  and  Eigenvalues  for  Seeded  Bands: 
Target  Area. 


variable 

EST  COHMUNaL 17 

T00I 

.93233 

7002 

.95521 

7  003 

.96|13 

TOO* 

.8692* 

T00S 

.90693 

TOO* 

.17188 

TOOT 

.72588 

T008 

.92377 

T009 

.96510 

T010 

.96113 

T011 

.96866 

7012 

.9617* 

7013 

.95583 

TOI* 

.9*592 

T01S 

.93763 

T016 

.9**32 

T01T 

.9*005 

T018 

.95105 

7019 

.95521 

TOZO 

.97*16 

T0Z1 

.9*016 

TOZZ 

.9*0*7 

T023 

.96866 

TOZ* 

.96352 

T02S 

.97*16 

T026 

•  90|  10 

T02T 

.93006 

T«2B 

.9526* 

7029 

.93763 

T030 

.93385 

T031 

.95583 

T  032 

.9526* 

T033 

.95372 

T  03* 

,89898 

T03S 

.97386 

T  036 

.9212* 

T037 

.97-156 

T038 

.9*8*3 

T039 

.9*67* 

T0*0 

.8*127 

TO*l 

.92372 

T0*2 

.9*67* 

T0*3 

.92)2* 

TO** 

.96073 

T0*5 

.96175 

T0*6 

.921*2 

T0*7 

.62323 

T0*8 

.957*2 

T0*9 

.9*8*3 

TOSO 

.97386 

T051 

,89579 

T052 

.92102 

T0S3 

.96082 

TOS* 

.89813 

TOSS 

.96o82 

T0S6 

.9*018 

T057 

.91003 

T  058 

.851*2 

T  059 

.9*917 

7080 

.93016 

T061 

•97| 06 

T  002 

.93302 

7003 

.9*888 

TOO* 

.93870 

TOOS 

.9**97 

7066 

.920  78 

7067 

.92811 

7068 

.9*888 

7069 

.88283 

70  70 

.89813 

TOTl 

.915*5 

7072 

.97108 

Component 


EIGENVALUE 

PCT  or  V AH 

CUM  PC 7 

53.85*30 

73.7 

73.7 

*.89*10 

6.8 

80.5 

A.Isoiv 

5.6 

86.2 

2.22267 

3.1 

89.3 

1.33385 

1.9 

91.2 

. 90319 

t.3 

92.4 

.06051 

1.2 

93.6 

.79256 

1.1 

9*. 7 

.1*292 

1.0 

95.8 

.58931 

.8 

96.6 

.51*57 

.7 

97.3 

.**061 

.6 

97.9 

.36978 

.5 

98.4 

.32229 

.* 

98.9 

.28169 

.* 

99.3 

.258*2 

.4 

99.6 

.23903 

.3 

100.0 

.21613 

.3 

100.3 

.19651 

.3 

100.5 

.17857 

.2 

100.8 

• i 52*6 

.2 

101.0 

.135*2 

.2 

101.2 

.13061 

.2 

101.* 

•  11363 

.2 

101.5 

.10813 

.2 

101.7 

.08765 

.1 

101.6 

.07853 

.1 

101.9 

.07502 

.1 

102.0 

.069*0 

.1 

102.1 

.06055 

.1 

102.2 

.05*98 

.1 

102.3 

.05367 

.1 

102.3 

.0*071 

.1 

102.* 

.03897 

.1 

102.* 

.03052 

.0 

102.5 

.028*3 

.0 

102.5 

.0269* 

.0 

102.6 

.02307 

.0 

102.6 

.0212* 

.0 

102.6 

.02071 

.0 

102.7 

.01382 

.0 

102.7 

.01222 

.0 

102.7 

.00930 

.0 

102.7 

.OOTil 

.0 

102.7 

.00530 

.0 

102.7 

.00396 

.0 

102.7 

. .00275 

.0 

102.7 

.0016* 

.0 

102.7 

-.00185 

•  •  0 

102.7 

-.003*7 

-.0 

102.7 

-.00**9 

-.0 

102.7 

—.00567 

-.0 

102.7 

-.00763 

-.0 

102.7 

-.00990 

-.0 

102.7 

-.01*02 

-.0 

102.7 

-.01817 

-.0 

102.6 

-.02002 

-.0 

102.6 

-.02353 

-.0 

102.6 

-.02613 

-.0 

102.5 

-.03130 

-.0 

102.5 

-.03759 

-.1 

102.5 

-.0*120 

-.1 

102.* 

-.05293 

-.1 

102.3 

-.06137 

-.1 

102.2 

-.08*79 

-.1 

102.1 

-.10293 

-.1 

ioz.o 

-.11227 

-.2 

101.8 

-.135*9 

-.2 

101.6 

-.iS7ii 

-.2 

101.4 

-.20420 

-.3 

101.1 

-.31952 

-.*  • 

100.7 

-.*98*6 

-.7 

100.0 

I 
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TABLE  A.3.C.  Coramunalities  and  Eigenvalues  for  All  Bands,  Non-seeded  Bands, 
and  Seeded  Bands:  Control  Area. 


ALL  BANDS 


variable 

EST  COMMUNAL  I 

cot 

.98636 

C42 

.9958s 

C03 

.99337 

CO* 

.97892 

COS 

.99533 

C06 

.99831 

COT 

.9969# 

COB 

.99828 

COO 

.96835 

CIO 

.99477 

Cll 

.99?64 

CI2 

.99572 

CIS 

.96346 

CIA 

.98381 

CIS 

.94|19 

CIO 

.99702 

C17 

.97763 

CIS 

.96544 

C19 

.97137 

Component 

EIGENVALUE 

PCT  Of  VAR 

CUM  PCT 

1 

14.452/5 

76.1 

76.1 

2 

1.2699# 

6.7 

82.6 

3 

.88598 

4.7 

87.4 

4 

.58713 

3.1 

90. S 

5 

•368VV 

1.9 

92.4 

6 

.34014 

1.8 

94.2 

7 

.2375/ 

1.3 

95.5 

6 

.193/1 

1.6 

96.5 

9 

.13059 

.7 

97.2 

10 

.107/4 

.6 

97.8 

11 

.89434 

.5 

98.3 

12 

.08635 

.5 

96.7 

13 

.07755 

.4 

99.1 

14 

.962/0 

.3 

99.4 

15 

.03569 

.2 

99.6 

16 

•02959 

.2 

99.6 

17 

.0234/ 

.1 

99.9 

18 

.01661 

.1 

100.0 

19 

.00032 

.0 

100.0 

NOT  SEEDED  BANDS 


VARIABLE  CST  COMMUNAL  I TY 


C01 

.79209 

C02 

.93852 

COS 

.89016 

CO* 

.89670 

COS 

.98033 

C06 

.97830 

C07 

.906*8 

COS 

.89256 

C09 

.85654 

CIO 

.68215 

Cll 

.84322 

C12 

.98033 

Cl  3 

.87393 

C14 

.67271 

CIS 

.78843 

C16 

.97830 

Cl  7 

.93852 

CIS 

.68188 

C19 

.89670 

Component 

EIGENVALUE 

PCT  Of  VAR 

CUM  PCT 

1 

14.4063V 

75.8 

7S.8 

2 

1 .40*68 

7.4 

83.2 

3 

.92517 

4.9 

88.1 

4 

•62445 

3.3 

91.4 

5 

.3848/ 

2.0 

93.4 

6 

.28338 

1.5 

94.9 

7 

.202V6 

1.1 

96.0 

8 

.iasol 

i.o- 

96.9 

9 

.13339 

.7 

97.6 

10 

• 10**6 

.6 

98.2 

11 

.09906 

.5 

98.7 

12 

.08615 

.5 

99.2 

13 

.07480 

.4 

99.6 

14 

•04712 

.2 

99.8 

!S 

.03867 

.2 

100.0 

16 

.01904 

.1 

100.1 

17 

.0128* 

.1 

ioo.2 

18 

•.00505 

-.0 

100.1 

19 

-.026// 

-.1 

100.0 

SEEDED  BANDS 


VARIABLE 

EST  COMMONALITY 

Component 

EIGENVALUE 

PCT  Of  VAR 

CUM  PCT 

COl 

.8*822 

1 

i4./iao6 

77.4 

77.4 

C02 

.91673 

2 

1.41624 

7.5 

84.9 

COS 

.69335 

3 

•8131V 

4.3 

89.2 

CO* 

.95527 

4 

.52588 

2.8 

91.9 

COS 

.97162 

5 

.36957 

1.’ 

93.9 

C06 

.98256 

.3400* 

|.8 

95.7 

C07 

.91655 

7 

.21433 

i.i 

96.8 

COS 

•9J53* 

8 

.20050 

1.1 

97 .6 

C09 

.90897 

9 

.12118 

•  0 

98.5 

CIO 

•95| 78 

10 
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TABLE  A.4.TA.  Station-Component  Correlationa  All  Bands:  Target  Area. 
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Table  A.4.TN.  Station-Component  Correlations  for  Non-Seeded  Bands: 
Target  Area. 
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Table  A.4.TS.  Station-Component  Correlations  for  Seeded  Bands: 

Target  Area. 
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Table  A.5.T.  Correlation  Matrix  for  the  Target  Area  over  All  Bands. 
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Table  A.5.T.  Correlation  Matrix  for  the  Target  Area  over  All  Bands. 
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TABLE  A.5.C.  Correlation  Matrix  for  the  Control  Area  over  All  Bands. 


®  N-  ®  ®  A  inn  lA  lAO  n  -N  O  «  CD  o 

-«  n  A  a  n  —  ©  \  A’——'*)—  — 

^  n  n* ©  n-  on  ^  -  nh  r>  —  m  0  —  t a 

m  O  AIM  •«  AJIM  M/)e  hi  -*«•  o  a  m  IMOO 

N-  ®  N-«  9>  OOO  CD  flOO  9  O0  N-  ®  ®  N-  CD  « 


N-  A  —  N*  —  ©  M  —  ©  A  M  —  0*  O  — «  A  O  — 

m  ©  —  o  ®  ct  ®  ®n*  m  ©  — 

O  O  O  O  A  O  ®  9  ®coMA®AA®A 

®  ®  nA  ®o®— oa®®— ooa  o  m  — 

®  r-  «  ®n*  ®  ^n.  ®  O®  ®  CD ®  a  lA^  ®  ® n- 


<n  no»N  ^  o®  o  —  a  ai  (Ao  -  on  ®  ®  n 

*  »  o»9  m«N  own  nn»  —  ®  ®  s-  —  o 

O  M  «—  r>  A®  o  <09  n  nn  o  ®  ®  CT*  MO 

OnOO®AOO  —  N-  ®  (AN.  CT*  —  (A  O  O  N. 

«  n-sn-n-ooooo®  o  coo  o  ia« 


maan-®n-o®mo®oo  O  ©O  CT  —  ® 

N-  9  <N  m  MO  N  OlA  o  OO  AO®  CM  O'  © 

®  n-  oo  n*  — ©  ®  on  o  nio«  tn  9- 

-9  0-0  — •©  >«M  N  OM  ®  A  —  M  c-n. 

NNNN*0>O  «  N.  ®  ®  N*  N*  O  |AO  **  N-® 


—  ®  O®  A  ON-  e  o A  ®  o®  o  AN.  M  Al- 

»N  N>  M®  —  OM  0  —  A  N-CT9  O  9®  A  -TO 

M  M  ON  O  ©  —  A  ♦  O  9  00  >*90  —  —  ®- 

A  O  99iS  O-  A  OM  n  09  A  O®  >  l/»© 

®  NN®NNO0«N«®NNN®0®0N 


®  ®  *»©  o  A®  ®  —A  «o«  9nn  Af'o 
®  n»  >»  —  o  —a  A  —  n  ct*  mo  ®  nn  ct  o ct 

o  OlN»<a  »-1in<*<9iAAO|OTOl/9' 
A  •  — M«lO  9®  OiO-'>-'Oi3  f®ni9N. 
A  N*  ®  N*N.  o  OO  ®<®0*N-  ON'NON’N  N-N. 


0  A  M©  ©  <x>  —  m  no  o  o®  o  no  n  an- 

O  A  MO  —  ®N>  ®  OCA  —  ® O  <9  O  —  O  o® 

o  ®  ®  o  (A  —  o  —  ® o  ®  on  >o  9®  n  o« 

«  A  OO  ®  ®  —  O  A  CM  r>  ®  CO  OCT*  A  A  O 

®  ®  N*©  N-  ON-  r»  N»®  N»  N-«  O  «  N-  ®  ®  CT 


o  n.  o—  ®  n—  o  — n-  mo©  ®  no  o 

—  AO  9  ©ON-  —  n.  ®  ©  O'—  A  O  O  O 

—  —  —  «  *  *  —  OAAOA®  —  ®  ®  O—  O 

CT  O  O  O  N.  o  0  0—0^0  <A  (M  CM  CM  O  A  O 

®N-®  CT  N-  N-  ®  ®  N>  ®  N»  —  N.  ®  ®  N-  ®  ®  o 


N-  O  —  A  O CM  —  ®  CT  ®  — O  CM  O  '  A  —  N-  O  <A 
—  0*®  0*0  O  —  O  —  A  *9  —  A*-  ®  AO  * 
— •  ®  A  ©  ®N»  ®M®—  OOAA*  ®  O  ©  — 
M  M  O  A  ®  A  —  O  M  ©  A  O  —  O  A  O  O  A 
®OA®N-®N-f-®®ANN®  ®  ®  A  O  ® 


f-®®  O  ACM  CT  ®  A  O  CM  N  A  CD'  -9  ®  ON-  O 

®  O®  CT  CT*  A  (MOM  O  —  ®  «*  CM1  O  AOA  <A 

O  OO  o  «  —  A  <A  A  ®  ®A  ®  O'A  CA  OCA  O 

O  <\J  ®  A  O  (A  MO-9  MOO  OO  M—  OO  O 

N  0*9  ®N>®  N.r»®N-®N'N>A  ANOA  ® 


>  A  O  CTON»aON>9—  ON-— M  AO®  —  ® 
000—0®  ®  9  N»  —  ®«M®AOA®A 
—  OO'O'OO  —  O  A  MOM  O'-  0*00©  ®  • 
ASM  CT*®  ®  —  AAMO®— 'O®©— A  M 
NNN  N-  N-CT*  O  9  N-  ®  ®N-  NS  A  ON-®  — 


O  M®  O  O  A  OO  —  AAON-CAOAOA  O 
M  ®  O  O  O  ®  ®  ®  ®  N.OM  O  O-  ©  A  O  — »  — 
—  *9  — *0*9  9  A  »i®i  —  M  9<9i9  O  >«A®*'  ©i 
Jt-O  O'O'O®  <  IV— *0- AM-O  O  *0  O  ®«MCT'M 

O  ®  0  9  ®  A  AAA  ®A®  ®  —  O  A  A®  9 


—  ®—  o  OO  ©— 00®©0O0M®©O 
®  M  O  O  ®0  A— 9  A— o  A  o  o  ®  MA  O 

n»  AO  ®  0O  O  O  A  OOO  OO®  N-  OA  — 

090  900  ®AO  O  O  M  900000  M 

A  A  ®  —  —  N-  ®  ®  A  N>  A  N.  ®  o  N-  N>  AO  ® 


N»  9'OM  N-  O A  O'— ®<0®®  —  ®-0  ®  -»0 
O*o*o  M  O'MO'O'O  A'O-N-O'O  OO  O  9  9* 
O  O  O®  O  O®  «*0000  0-0  O  ®  — 
—  O  O  O  M  ®  O  00000-0  O  O  M  ®  O  O 
o«ooo®oo«o®oo  ®  ® O  «  A  ® 


OOOAO®AOAOOAO®MA®0— 
o  OOAOO®  CT*  —  A  AMA  M  «  O  O  CT*  A 
0*0  N-OlOiNN-M  >O-A9O*A  O0  Ol»  — * 
7*0  —  A  —  ®®l0'®0  —  9— *9-0  9  M  M O 
O  O  ®0  «  OO  ®  9®  O  ®0  A  ®0  >  ®0 


OOO®  ®  •— (M  OO®  0®A  OOO  o  o  — 
O  O  OO  o  OO  ®  ®  —  —  ®  A  ®  M  O  ®  — — 
®  O  ®e  O  M®  o  OO  CM  OO  o  O  —  O  — — 
OO— ®AA  —  OO  -  —  MA  OAAOMO 


AO®  ®  99  O-OO-n-  ®  O  O  t  O  — AM  M 
AAO’OO®  0  0  0*0  OO  OAO'MO—  O 
OOO  O  A  O  OOMOO®  09®  O  « A  O 
A  O  O  OOO  MO««««OM  0900®  —  o 
®  OO  00009®  ®o®  o®®oooo 


®  A®  OO®  «  A  —  —  ooooooooo 
®MO  «  MO  ®  O  ®  o  M  O  OOMOO®  ® 
Ot®>>  AO  M09  O  A  9  9  9  OlM  AO*A 
MOOiO  — ®i>*AO>  —  — O  M  M  0*91®®  —  • 
0®0  OOO  OO®  OOO  ®o®oooo 


ortoo®®  ®CVCMOOO®®A®M  —  — 

—  AO  —  OO  OOO  MOM  0—0  9  —  A  o 
M  AO  O  O®  OO®  —  OAOOMO  ®0  • 

—  — «  OO  —  0®®MO  —  0®M0«A  O 
00®00®  9®®®OOOAA®®AO 


—  M  n  ®  A  ®o  OOO  —  Mn  O  A  ®  o  ®  o 
O  O  OO  O  OO  o  o—  —  —  —  —  —  —  —  — — 

u  uuu uuu o  ou u  uu u  uu uuu 


—  Mn  OA®  0«0  o— MnoA®0«  o 
OOO  OOO  OOO  —  —  —  — —  —  — —  — 

o  UO  U  OU  UJO  u  uo  ucuuuuu 


Report  Numbers 
M388,  ONR-110 

M409,  ONR-111 
M410,  ONR-112 
M417,  ONR-117 

M419,  ONR-118 
M420,  ONR-119 
M428,  ONR-122 
M440,  ONR-126 

M442,  ONR-127 
M467,  ONR-133 

M469,  ONR-134 

M490,  0NR-13S 

MS 14,  ONR-147 
MS21,  ONR -149 


1 


OTHER  TECHNICAL  REPORTS  ON  THIS  CONTRACT 


Hanson,  Morgan,  Bach,  Charles  L.  and  Cooley,  Edward  A., 
Bibliography  of  Statistical  and  Meteorological  Methodology 
in  Weather  Modification,  September,  1976. 

Bradley,  Ralph  A.  §  Srivastava,  S.S.,  Correlation  In 
Polynomial  Regression,  March,  1977. 

Bach,  Charles  L.,  An  Interpretive  History  of  30-years 
(1945-1975)  of  Weather  Modification,  March,  1977. 

Bradley,  Ralph  A.,  Srivastava,  S.S.  and  Lanzdorf,  Adolf, 
Summarization  of  Precipitation  Data  In  a  Weather  Modification 
Experiment:  I.  A  Response  Surface  Approach,  June,  1977. 

Gleeson,  T.A. ,  Data  Summarization  In  a  Weather  Modification 
Experiment:  II.  Concomitant  Variables,  June,  1977. 

Hanson,  Morgan  A.,  Rank  Tests  in  Weather  Modification  Experi¬ 
ments,  June,  1977. 

Serfling,  R.  J. ,  Toward  a  Nonparametric  Covariance  Analysis 
of  a  Weather  Modification  Experiment,  August,  1977. 

Hanson,  Morgan  A.,  Barker,  Lawrence  E. ,  and  Hunter,  Charles  H. 
Bibliography  of  Statistical  and  Meteorological  Methodology  in 
Weather  Modification,  II,  October,  1977. 

Scott,  Elton,  Data  Summarization  in  a  Weather  Modification 
Experiment:  III.  A  Multivariate  Analysis,  June,  1978. 

Bradley,  Ralph  A.,  Srivastava,  Sushil  S.,  and  Lanzdorf,  Adolf, 
An  Examination  of  the  Effects  of  Cloud  Seeding  in  Phase  I  of 
the  Santa  Barbara  Convective  Band  Seeding  Test  Program,  June, 
1978. 

Hanson,  Morgan  A.,  Barker,  Lawrence  E.,  Bach,  Charles  L. , 

Cooley,  Edward  A. ,  and  Hunter,  Charles  H. ,  A  Bibliography 
of  Weather  Modification  Experiments,  July,  1978. 

Bradley,  Ralph  A. ,  Srivastava,  Sushil  S.  and  Lanzdorf,  Adolf, 
Some  Approaches  to  Statistical  Analysis  of  a  Weather  Modifi¬ 
cation  Experiment,  January,  1979. 

Scott,  Elton,  A  Multivariate  Methodology  for  the  Analysis  of 
Weather  Modification  Experiments,  August,  1979. 

Bradley,  Ralph  A.,  Scott,  Elton,  Randomization  Tests  in  Support 
of  Some  Statistical  Analyses  of  a  Weather  Modification  Experiment, 
October,  1979. 


67 


UNCLASSIFIED _ _ _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 

1.  REPORT  NUMBER  2.  GOVT.  ACCESSION  NO. 

FSU  Report  No.  M514 

ONR  Report  No.  147 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  {and  subtitle) 

A  MULTIVARIATE  METHODOLOGY  FOR  THE  ANALYSIS 

OF  WEATHER  MODIFICATION  EXPERIMENTS. 

5. ‘TYPE  OF  REPORT  6  PERIOD  COVERS! 

Technical  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

FSU  Report  M514 

7.  AUTHOR  (a!) 

Elton  Scott 

8.  CONTRACT  OR  GRANT  NUMBER^s) 

ONR  N00014-76-C-0394  ^ 

~T.  PERFORMING  ORGANIZATION  NAME  ANl)  AbDRESS 

The  Florida  State  University 

Department  of  Statistics  ' 

Tallahassee,  Florida  32306 

10.  PROGRAM  ELEMENT,  PROJECT,  TASS 
AREA  &  WORK  UNIT  NUMBERS 

HI  CONTROLLING  0FPl6E  NAME  ANb  ADDRESS 

Statistics  and  Probability  Program 

Office  of  Naval  Research 

Arlington,  Virginia  22217 

12.  REPORT  DATE” 

August,  1979 

13.  NUMBER  OF  PAGES 

67 

, 

14.  MONITORING  AGENCY  NAME  &  ADDRESS (If  different 
from  Controlling  Office) 

15.  SECURITY  CLASS (of  this  report) 

Unclassified 

15a. declassification/downgrading 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT (of  this  Report)  - 

Approved  for  public  release;  distribution  unlimited. 

17.  DlSlklbbilun  siAitrusra  lot  the  abstract  entered  in  Block  20,  if  different  from 
report) 

1ST  SUPPLEMENTARY  NOTES  " - — - 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Weather  modification,  cloud  seeding,  principal  components,  summarization  of  data 
dimension  reduction,  Hotelling's  T^  statistic,  cloud-seeding  effects,  data  analysis 
applied  regression,  multivariate  analysis.  ’ 

20.  ABSTRACT  (Continue  on  reverse  elds  if  necessary  and  identify  by  block  number} 

UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE _ _ 

REPORT  DOCUMENTATION  PAGE 

20.  ABSTRACT 

This  paper  develops  applications  of  multivariate  statistical  models,  particu¬ 
larly  principal  component  analysis,  to  the  analysis  of  data  from  weather  modification 
experiments.  The  efficacy  of  these  multivariate  applications  is  examined  by  apply- 

A 

ing  the  proposed  models  to  data  from  Phase  I  (1967-71)  of  the  Santa  Barbara  Convective 
Band  Seeding  Program  conducted  for  the  Navy  by  North  American  Weather  Consultants. 
Multivariate  summary  measures  of  precipitation  are  developed  and  multivariate  methods 
are  given  to  analyze  the  effects  of  cloud-seeding  on  precipitation.  Results  from 
these  models,  based  on  the  above-mentioned  data  set,  are  reported  along  with  conclu¬ 
sions  and  suggestions  for  further  work.  An  appendix  provides  detailed  summary  sta¬ 
tistics  for  the  analyses. 

I 


